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Abstract 

Entropy production in the compression stage of heavy ion coUisions is discussed 
within three distinct macroscopic models (i.e. generalized RHTA, geometrical overlap 
model and three- fluid hydrodynamics). We find that within these models ~80% or 
more of the experimentally observed final-state entropy is created in the early stage. 
It is thus likely followed by a nearly isentropic expansion. We employ an equation 
of state with a first-order phase transition. For low net baryon density, the entropy 
density exhibits a jump at the phase boundary. However, the excitation function of 
the specific entropy per net baryon, S/A^ does not reflect this jump. This is due to the 
fact that for final states (of the compression) in the mixed phase, the baryon density 
Pb increases with y/s, but not the temperature T. Calculations within the three- fluid 
model show that a large fraction of the entropy is produced by nuclear Shockwaves in 
the projectile and target. With increasing beam energy, this fraction of S/A decreases. 
At y/s = 20 AGeV it is on the order of the entropy of the newly produced particles 
around midrapidity. Hadron ratios are calculated for the entropy values produced 
initially at beam energies from 2 to 200 AGeV. 



Work supported bx) BMBF, DFG, GSI aud the German Academic Exchange 
Sermce (DAAD). 
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1 Introduction 



A large amount of entropy can be produced in the initial off-equilibrium stage of energetic 
nuclear collisions [||. The subsequent expansion is, on the other hand, often assumed to 
be isentropic. The linear increase of the pion multiplicity with beam energy observed at 
the BEVALAC is consistent with a picture 0, where more than 80% of the specific 
entropy is produced initially. In this scenario the entropy produced during the compression 
stage is closely linked to the finally observable relative particle yields. 

In the present paper, the entropy produced within the compression stage of highly rel- 
ativistic heavy ion collisions is calculated within three distinct macroscopic models (i.e. 
generalized RHTA, geometrical overlap model, three- fluid hydrodynamics). Several compar- 
isons of experimental data on relative hadron multiplicities with calculations of chemically 
equilibrated hadronic gases indicate that the entropy per net baryon at (chemical) freeze-out 
is on the order of S/A = 10 ± 2.5 at AGS and S/A = 40 ± 10 at SPS energies §, 0, g ||. 
The microscopic transport model UrQMD also exhibits a chemical composition in the cen- 
tral volume of Au+Au collisions at AGS energies which is compatible with an ideal hadron 



gas with S/A= 12 ||T0|. Moreover, S/A values in this range yield the measured rapidity 
distributions and transverse momentum spectra of produced hadrons at AGS and SPS, if an 
isentropic expansion of the hot and dense matter created during the early compression stage 
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The question arises, whether and how such a tremendous amount of entropy can indeed 
be created in the initial stage of the reaction, before local equilibrium is established. If true, 
this would speak in favour of the picture that the compression stage produces most of the 
finally observed entropy and that it is followed by a nearly isentropic expansion. 

The present study focuses on the excitation function of S/A, i.e. its dependence on the 
bombarding energy. In particular, it is interesting to study the behaviour of S/A for energies 
where the matter in the central region enters the mixed phase and the pure QGP phase. 
This analysis reveals whether the produced entropy exhibits an "anomalous" behaviour which 



could serve as a direct sign of the QGP phase transition |T2|. 

Also, the fraction of entropy produced both by nuclear Shockwaves in the projectile 
and target, on the one hand, and by particle production (and their thermalization) around 
midrapidity, on the other hand, is studied as a function of bombarding energy within the 
three-fluid model. 



2 Models for the compression stage 

In this section, the equation of state and the three collective models which will be employed 
in section ^ to calculate entropy production are described. 

2.1 The Equation of State 

The equation of state (EoS) of the highly excited nuclear matter must be specified in order to 
determine entropy production in the compression stage. Here, the a—uo model is used for 
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the interactions between the nucleons. The parameters of the Lagrangian are fitted such that 
the properties of infinite nuclear matter in the ground state are reproduced. The pressure of 
an ideal pion gas in chemical equilibrium is added to the pressure of the interacting nucleons 
to account for thermal pion production. Other particles like A's and p's are not included 
in the EoS, since they probably do not reach thermal and chemical equilibrium in the early 
stage of the reaction. 

At high temperatures or baryon-chemical potentials, a phase transition to the quark- 
gluon plasma (QGP) is assumed to occur. Here, the MIT bag model is used to calculate the 
EoS fl^ for massless, non-interacting gluons and u, d quarks. The pressure as a function of 
temperature T and quark-chemical potential Hq reads [[Tsl] 
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A bag parameter of B^/'^ = 235 MeV is used. The two EoS are matched by Gibbs' conditions 
of phase equilibrium. The resulting EoS thus exhibits - by construction - a first-order phase 
transition at a critical temperature (for vanishing pb) of Tq ~ 170 MeV. Baryon-, entropy-, 
and energy-densities can be calculated via the usual thermodynamic relations 



PB 



dp 



(2) 



where the baryochemical potential 
refer to 



d/iB 
dp 
df 

IJ-B 

Ts -p + hbPb , 

= 3/ig. For a more detailed discussion of this EoS 



For this EoS, the entropy per net baryon s/pb exhibits a jump at the phase transition 
temperature, if pb is not too large (Fig. |I]). However, the excitation function of s/pb does 
not show a corresponding jump at a given specific collision energy, as will be discussed in 
section RTT. 



2.2 Geometrical Overlap model 

The first model discussed is a simple schematic geometrical overlap model, where the two 
Lorentz-contracted nuclei are placed on top of each other. Then the following baryon density 
Pb and energy density e are obtained: 

Pb = 2po7cM (3) 
e = VspolCM = ^pI ■ (4) 

Note that the compression stage is not treated dynamically within this model. Therefore 
the compression achieved is independent of the EoS. However, the EoS determines the tem- 
perature and the entropy density of the compressed matter. 
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Figure 1: Entropy per net baryon as a function of temperature at various fixed net baryon 
densities for the EoS used here (see text). 



2.3 One-fluid hydrodynamics 

We now turn to the one-fluid hydrodynamical model, restricting the calculations to an ideal 
fluid, i.e. viscosity and thermoconductivity effects are neglected [|, |17]. 



In this model, compression and entropy production occur due to shock discontinuities, 
or by a combination of shock- and isentropic compressional waves if the compressed mat- 
ter is thermodynamically anomalous (as e.g. matter in the mixed phase) [|TB], |T^, 
The assumption of instantaneous thermalization leads to the maximum energy- and baryon 
number-deposition in the central region that is consistent with energy-, momentum-, and 
baryon number-conservation. 

As an example. Fig. ^ shows a three-dimensional one-fluid dynamical calculation of the 
collision of two Au-nuclei at a/s = 4.9 AGeV and vanishing impact parameter, 6 = 0. The 
upper panels of Fig. || show reaction plane contour plots of the baryon density 7Pb/po in 
the center of mass frame. Two Shockwaves can be clearly observed as high density gradi- 
ents. They develop in the center and propagate outwards. The entropy produced by these 
Shockwaves is depicted in the lower panel. The entropy per net baryon is computed as 

A J d-^x'jpB 

where s is the entropy density (cf. eq. The compression stage is finished once the shock 
fronts reach the back sides of the Lorentz-contracted nuclei. Then an isentropic expansion 
follows (remember that we assume a non- viscous fluid). 

The specific entropy (and also the energy density, the pressure, etc.) of the central region 
(close to x,y,z = 0) in relativistic nucleus-nucleus collisions at 6 = can also be computed di- 



rectly by solving the Rankine-Hugoniot-Taub Adiabate (RHTA) equation |T^, [T^, |T 
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Figure 2: Time evolution of baryon density and S/A in a 3D one-fluid hydrodynamical 
calculation of Au-^Au a,t ^/s ^ 4.9 AGeV. 
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thus assuming a one- dimensional "slab-on-slab" collision. The RHTA determines for a given 
initial state and EoS all shock wave solutions consistent with energy and momentum con- 
servation. Matter in the mixed phase, however, is thermodynamically anomalous and single 
shock waves become mechanically unstable. Therefore, construction of the generalized RHTA 



is required |T^, |T^ . In the generalized RHTA, the unstable shock wave solutions are replaced 
by a combination of shock waves and an isentropic compression wave. For matter in the pure 
QGP phase, a single shock becomes again stable. For a given bombarding energy, the gen- 
eralized RHTA thus determines the thermodynamic properties of the compressed matter 
"final" state). S/A is independent of i/i for all bombarding energies with "final" states in 



the phase coexistence region (cf. section SJ^), due to the fact that a simple compression wave 
does not produce any additional entropy. 

2.4 Three-fluid hydrodynamics 

The assumption of instantaneous local thermalization in the one-fluid model prevents the 
expected partial interpenetration of projectile and target. A more realistic description, 
allowing for finite mean free path effects, is obtained if two distinct fluids for projectile and 
target are introduced, which interact locally [p2| , These interactions lead to energy and 



momentum exchange and to a smooth deceleration of the projectile and target fluids. At 
AGS and SPS energies, it is reasonable to introduce a third fluid that "collects" the produced 



particles around midrapidity ||2^, The model assumes that the fluids are well separated 
in rapidity (but, of course, not in coordinate space) during the early stage of the reaction. 
Binary particle collisions increase the overlap of the fluids in momentum space and thus drive 
the system into (local) mutual kinetic equilibrium. At this point, the various fluids merge 
into a single fluid [^, ^ and the subsequent expansion proceeds as in the one-fluid 
hydrodynamical model. The time scale for this kinetic equilibration of the central region 



was estimated in refs. ||2J, ^ to be ~ 2R/^cm- Here the interactions between projectile 
and target nucleons have been regarded as a sum over incoherent binary nucleon-nucleon 
scatterings (with cross sections as in vacuum). Since we focus on the early (compressional) 
stage of the reaction, the one-fluid transition is not discussed here. 

The most essential difference to the one-fluid model is that the shock fronts are smeared 
out considerably, i.e. projectile and target interpenetrate strongly. Consequently, the baryons 
are less compressed [^, ^ in the three-fluid model. Most of the energy loss of the incoming 
nucleons results in particle production at midrapidity, i.e. a large part of the energy loss 
of the nucleons is transfered to the mid-rapidity regime (third fluid) and only a smaller 
fraction of the energy is deposited in the baryon fluids. The physical mechanism for entropy 
production is very different in the three-fluid model than in the one-fluid hydrodynamical 
model, therefore the ratio of the thermal to the compression energy at midrapidity differs 
from that given by the RHTA equation. 
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3 Results 



We now turn to the discussion of entropy production as a function of the bombarding energy 
for the three different models. 

3.1 Excitation function of Specific Entropy 

The excitation function of the entropy per (net) baryon is shown in Fig. ^ for the three 
different models. 

1. The geometrical overlap model assumes that energy and baryon density increase as '^'qm 
and 7cM, respectively (cf. eqs. (0^)). As a consequence, S/A increases continuously 
with energy. The compressed matter enters the mixed phase at y/s k, 2A AGeV with 
a specific entropy of S/A = 4.5 . Pure QGP matter is produced with S/A = 12 at 
^/s ^4.6 AGeV. However, while sweeping through the mixed phase by increasing the 
bombarding energy, pB and S/A increase smoothly, while T is roughly constant, cf. 
Fig. ^. Thus, the entropy jump seen in Fig. |1] is not reflected in the S/A{^/s) excitation 
function. 
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Figure 3: Excitation function of entropy per participating net baryon as calculated within 
the various models. For the three-fluid model, both the entropy of projectile- and target 
fluids only (1+2), as well as that of all three fluids (1+2+3) are shown. 



2. In contrast, in the shock model (ID RHTA) S/A first increases with the beam energy 
and reaches a plateau at ^/s ^2.7 AGeV. At this beam energy the matter enters the 
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mixed phase. Single shock fronts cease to be the stable solution []TB|, |TB|, jTU], ^ in 
this thermodynamically anomalous region. Rather, further compression can only be 
achieved by a simple compressional wave which, however, conserves the specific entropy. 
Thus, even when the bombarding energy is increased, S/A remains constant (isentropic 
compression) until the anomalous region is left again. Hence, the baryon density 
increases with energy while the temperature decreases, as seen in Fig. H. In this model. 




Figure 4: Temperature and (net) baryon density of the compressed matter within the overlap 
model and the shock model (generalized ID RHTA), respectively. The phase coexistence 
region is also indicated. 



the participant matter reaches the pure QGP phase at a temperature of T = 65 MeV 
at ^/s ~ 3.3 AGeV, which is significantly less than in the overlap model (Fig. The 
specific entropy then increases monotonically with energy. In the AGS energy range, 
where a single shock is again stable, we find S/A 10, in agreement with the numbers 
obtained from the data analysis on hadron yield ratios [^]. However, the energy density, 
baryon density, and entropy density are e = 5 GeV/fm'^, pB = 2.0 fm~^, s = 19 fm~^. 
These values seem unreasonably high at this energy. These huge densities and the very 
early phase transition to pure QGP matter are due to the assumption of instantaneous 
stopping and thermalization p2|, |2^, At higher energies, the discrepancies become 
even more pronounced: e.g. for Pb+Pb at SPS, e ~ 50 GeV/fm^, ^ Q fm~^, and 
s ~ 150 fm~^. These values are close to the analytic expressions e/eo = 47^j^^ — 3, 
pg/po = A-fcM - ^hcM (eo = Po ■ 923 MeV= 148 MeV/fm^, po = 0.16 fm'^) derived 
from the RHTA with an ultrarelativistic EoS, p = e/3. Because these densities are so 
huge, the specific entropy is rather low, S/A = 25. This is in fact significantly below 
the value 40 ± 10 discussed in [|], |^ H, ^ |1T|. Slightly larger S/A values are found 
when the number of degrees of freedom in the QGP is increased, e.g. by including 
strange quarks. However, at the SPS the compression phase in this model lasts only 
for 1.1 fm/c (the time when the Shockwaves reach the back of the nuclei). Within 
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this very short time, strange quarks may not reach thermal and chemical equilibrium. 
In the subsequent expansion, strange quarks (resp. strange hadrons) should be taken 
into account. 
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Figure 5: Time evolution of the entropy per net participant baryon for Ph + Ph reactions at 
various energies, calculated within the three-fluid model. 



The three-fluid model predicts the entropy per net participating baryon as a function 
of CM-time as depicted in Fig. |^. In contrast to our earlier work [^, ^ 27|, in this 



calculation the EoS described in section is employed for all three fluids. The third 
fluid is net baryon free since we assume that the direct baryon transfer to midrapidity 
is small in the flrst collisions with large rapidity gap [^. To calculate SjA within the 
three-fluid model a temperature cut of T > 50 MeV has been applied to determine 
participant baryons. This corresponds to a mean transverse momentum of (pt) ~ 270 



MeV, roughly in accordance with the cut employed in [jl2| . Total entropy and baryon 
number have then been calculated by summing over all three fluids (resp. projectile 
and target fluids only for baryon number). Thus: 



S 



E /d3a;7se(T-50 MeV) 

fluids 1,2,3 /„\ 

(6) 



fluids 1,2 

At early times, t < 2R/'ycM, kinetic equilibrium between the fluids is not established. 
Thus, summing over the entropy and baryon currents of the individual fluids as in 
eq. seems physically more meaningful than calculating the entropy and baryon 
densities from the total energy-momentum tensor T^'^ = T/"^ + Tg^'' + T^" and total 
baryon current = jf -|- assuming that they are of the form 

T^"" = {e + p)uV -pgf"" , (7) 
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f = PBU^ 



as appropriate for a single (ideal) fluid [|T^]. Fig. || shows that S/ A saturates rapidly. 



and is essentially time independent for later times. This also speaks for the fact that the 
numerical "production" of entropy is negligible in the present calculations, which focus 
on the early compression phase of nuclear collisions. Baryon number A is conserved 
with an accuracy of ~ 0.1% within our model, and the time independency of S/A 
shows that this holds also true for the entropy 5, at least for the early stages regarded 
here. 

This plateau value is plotted in Fig. | as a function of ^/s (diamonds). If one omits 
the third fluid in calculating the entropy, S/A saturates similarly. The plateau values 
reached in this calculation are also depicted in Fig. ^ (squares). One observes that 
with increasing bombarding energy a larger fraction of the total entropy results from 
thermalization of the particles produced around mid-rapidity. At SPS, 50% of the 
total entropy is due to Shockwaves in the projectile and target fluids. After ~ 2 fm/c 
(measured in the CMS) the central region is in equilibrium and the isentropic 
one-fluid expansion sets in with S'/A = 38 (for y/s = 20 AGeV) and S/A = 35 (for 
y/s = 18 AGeV), respectively. This value is 10 units above the one-fluid result quoted 
above. Consequently, particle ratios involving antibaryons (like N/tt, N/N, A/A) are 
expected to increase considerably, while meson-meson and baryon-baryon ratios (like 
K/tt, A/N) do not change very much ^6 . 



Fig. ^ depicts the fluid rapidity density of the entropy within the various fluids. The 
fluid rapidity is defined as 

tanhi] = v\\ , (9) 

where v\\ is the fluid velocity component parallel to the beam axis. One observes 
that the largest contribution to the entropy around ?7 = is due to production of a net 
baryon free QGP within the third fluid. The entropy produced by the Shockwaves in the 
projectile and target is found at i] = ±1. Note, however, that the entropy distribution 
in the final state will be different from that depicted in Fig. ^. The high pressure 
gradients lead to reacceleration of the strongly decelerated fluids, i.e. expansion sets 
in. As shown in , the resulting rapidity distribution of the pions, including thermal 
smearing of the particle momenta, is Gaussian-like and does not exhibit peaks (like 
dS/dr] at tcM = 3 fm/c, cf. right panel of Fig. |^). In other words, at SPS energy the 
three-fluid model does not predict boostinvariant initial conditions for the expansion 
stage. Thus, neither the further evolution of the various rapidity bins nor their final 
break-up into hadrons (cf. section |3^) proceeds independently. This is in contrast to 
the dynamics at collider energies {^/s > 200 AGeV) as discussed in 



3.2 Particle Ratios 

The initial compression stage is followed by a nearly isentropic expansion, until the freeze- 
out is reached. This scenario allows to calculate particle ratios as follows. We assume that 
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at (chemical) freeze-out the net baryon density is p^^ = pq/2 and the net strangeness of the 
system is zero. The third thermodynamic variable, the specific entropy, is calculated within 
the three-fluid model, cf. section These quantities determine the chemical composition 
of the fireball - in local thermodynamical equilibrium - unambiguously. 

Investigations of hadron ratios at AGS and SPS energies indicate that not only pions 
and nucleons, but also heavier hadrons (including strange particles), and even clusters like 
deuterons are close to chemical equilibrium @, ^ ||, |^ at freeze-out. Therefore an ideal 
hadron gas model with all known hadrons and resonances up to 2 GeV/c^ mass is employed 
to determine the observable hadron ratios . Feeding from post freeze-out decays of heavy 
resonances is also taken into account. 

We thus imphcitly assume that the (fast) compression stage is followed by a rather long 
expansion stage before (chemical) freeze-out occurs. If this is indeed the case, it is reasonable 
that the EoS valid within the first fm/c is different from that at the late freeze-out stage. 
In principle, the chemical equilibration of the heavier hadrons and resonances can produce 
additional entropy, which is however neglected in the present analysis. 

It has been pointed out that the observed number of pions per net baryon riT^/uB exhibits 
an interesting behaviour |12|: For collisions of nuclei in the AGS energy region {^/s = 5 — 6 
AGeV), the ratio n-j^/nB is smaller than that in p + p reactions at the same energy (per 
nucleon). At SPS energies (a/s = 18 — 20 AGeV), however, the difference of the ratios is 
positive. Fig. |^ shows that this can be attributed to the much higher specific entropy at the 
SPS as compared to the AGS energy. Let us point out again that the excitation function of 
this difference changes smoothly - it does not exhibit a jump at some specific bombarding 
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Figure 7: The number of pions per net participant baryon in Pb+Pb minus that ratio for 
p+p as calculated from the S/A values obtained from the three-fluid model. The feeding 



due to decays of resonances is taken into account. Experimental data [|T^] are also shown. 
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energy (cf. sections ^.l] , as one might expect from Fig. |I[ The average number of pions 
in p + p reactions is parametrized as a function of as ||30| 



= -1.7 + 0.841ns + l.Os"^ (10) 



1 



= -2.6 + 0.871ns + 2.7s~2 (11) 
3 

= -{n^+ + n^-) . (12) 

At ^/s ~ 20 AGeV, about 50% of S/A are due to thermalization of the energy loss in 
the third fluid around mid-rapidity (at SPS, this fluid is in fact in the QGP phase with 
fig = fis = 0) and about 50% of S/A are due to nuclear Shockwaves in the projectile and 
target fluids (which are baryon-rich, but are nevertheless also in the QGP phase). Both 



contributions therefore are essential to understand the data. As discussed in section 
the three fluids, however, reach kinetic equilibrium later on and merge into a single fluid. 
Therefore we consider a single fireball as the source of the frozen-out hadrons, rather than 
letting each of the fluids break up independently. 

The hadron ratios at AGS and SPS, cf. Fig. ||, are quite close to the data discussed 
in the literature f^, |^, ||, ||, |10|, For such a simpleB estimate of hadron production in 



nuclear collisions, deviations from the experimental ratios by up to factors of two have to be 
expected. Nevertheless, it is clear from Fig. |^ that the simultaneous measurement of various 
hadron ratios, like vr/ {j3 — B^, d/N and, in particular, B/B (provided antibaryons also 
reach chemical equilibrium) allows to determine the produced entropy in the energy range 
between the AGS and the SPS. In contrast, the K/n-iaiio is practically constant. 



4 Summary and outlook 

In the present paper we have calculated the amount of entropy produced in the compression 
stage of relativistic heavy ion collisions using three different macroscopic models. 

The three-fluid model is the only of these models which accounts for (kinetic) non- 
equilibrium effects between the mutually interacting fluids. Calculations within this model 
show that the amount of entropy which is due to production of a third fluid (in addition to the 
projectile and target fluids) of thermalized secondaries around midrapidity can - at higher 
bombarding energies - exceed the amount of entropy produced by nuclear Shockwaves in 
the projectile and target nuclei. The total specific entropy S/A produced within this model 
is consistent with the S/A values extracted from data using relative particle yields from 
equilibrated hadron gases. We find, e.g., S/ A = 11 for AGS and S/ A = 38 for SPS energies. 

^We neglect e.g. that some hadron species might not reach their chemical equilibrium abundance or 
decouple earlier (at higher than others. We also do not discuss here the rapidity dependence of hadron 
ratios (cf. e.g. for a calculation within the microscopic transport model UrQMD) to study whether 
the various contributions to dS'/d?], as depicted in Fig. ^, can be disentangled. This requires the one-fluid 
solution in that part of the forward light-cone between the freeze-out hypersurface and the hypersurface 
where local kinetic equilibrium between the fluids is established. Such a calculation is out of the scope of 
the present work. 
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Figure 8: The excitation function of various particle ratios as calculated from the S/A 
values obtained from the three-fluid model. Feeding due to decays of resonances is taken 
into account. 
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The present calculations therefore support the picture that most of the entropy production 
in relativistic heavy ion collisions occurs in the early compression stage. Then a (nearly) 
isentropic expansion follows. 

The excitation function of the specific entropy S/A[y/s) docs not exhibit any threshold 
signatures of the phase transition to the QGP, which is included in the EoS. Nevertheless, 
the phase transition might still be visible in the excitation functions of other observables, 
like e.g. directed in-plane flow. 

The predicted particle yields exhibit a smooth excitation function, too. Simultaneous 
measurement of various hadron ratios allows to determine the amount of entropy produced 
in relativistic heavy ion collisions. 

Future work will also consider the rapidity and transverse momentum spectra of various 
particles at freeze-out. A different EoS will be used during the expansion stage in order to 
account for the chemical equilibration of additional particle species. The evolution of the 
rapidity density of the entropy until freeze-out will be studied. 
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